Low-lying energy spectrum of the cerium dimer 
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The electronic structure ol Ce2 is studied in a valence bond model with two 4/ electrons localized 
at two cerium sites. It is shown that the low-lying energy spectrum of the simplest cerium chemical 
bond is determined by peculiarities of the occupied 4/ states. The model allows for an analytical 
solution which is discussed along with the numerical analysis. The energy spectrum is a result 
of the interplay between the 4/ valence bond exchange, the 4/ Coulomb repulsion and the spin- 
orbit coupling. The calculated ground state is the even = A = E = level, the lowest excitations 
situated at ~ 30 K are the odd ^ = A = E = state and the ^65 doublet = ±5, A = ±6, E = =f1)- 
The calculated magnetic susceptibility displays different behavior at high and low temperatures. In 
the absence of the spin-orbit coupling the ground state is the ^E~ triplet. The results are compared 
with other many-electron calculations and experimental data. 

PACS numbers: 31.15.xw, 31.15.vn, 33.15.Kr 
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I. INTRODUCTION 

Understanding the electronic nature of chemical bonds 
of /-elements is an issue which has been attracting wide 
attention of investigators working in diverse areas of re- 
search. The main question is the involvement of 4/ 
states in the formation of chemical bonds in molecules 
and in solids. In general, one can speak of two ex- 
treme scenarios. That is, either /—electrons partici- 
pate in chemical bonding which is the case for 5/-states 
in the U2 molecule [l|, or they stay essentially local- 
ized at each site forming core-like states like in lan- 
thanides. This phenomenon known as the localization- 
delocalization paradigm in solid-state physics [2!| is the 
subject of intense theoretical considerations and discus- 
sion. 

In various compounds cerium is believed to display 
both types of 4/ behavior which makes it the key ele- 
ment in these studies. Indeed, in solids Ce and its com- 
pounds exhibit a broad range of unusual electronic, mag- 
netic and structural properties called correlation effects 
3. While much effort has been invested in understand- 
ing of the complex electronic behavior in sophisticated 
cerium compounds, very little attention has been given 
to the simplest cerium chemical bond 04^1 j which is re- 
alized in Ce2. Meanwhile, owing to the progress of the 
spectroscopy of small metal clusters Q precise experi- 
mental data on the cerium dimer Q have been reported 
recently, followed by first reliable many-electron calcula- 
tions 

Cao and Dolg Q have performed a series of cal- 
culations of Ce2 using the complete active space self- 
consistent field (CASSCF) method, multireference con- 
figuration interaction (MRCI) method and coupled clus- 
ter approach [CCSD(T)] and concluded that its ground 
state "superconfiguration" is 4/^4/^(6sCTg)^(5(i7r„)'*. Six 
low- lying molecular terms ^S~, '^Sg , ^6g 

and ^6u) are "virtually degenerate" so that the exact 
refinement of the ground state is very difficult. The au- 



thors then restricted themselves to the and "^S^ 

ground states for which in the scalar relativistic approx- 
imation they obtained the following spectroscopic con- 
stants: Re = 2.62 ± 0.02 A, De = 1.60 ± 0.41 eV and 
ujc = 201 ± 13 cm~^. The low-lying excitations were at- 
tributed to the core-like properties of 4/ states. In the 
spin-orbit calculation the authors obtained the [ft = 0]'^^ 
ground state followed by two low- lying levels 0^ (< 1 
cm~^) and 5„ (< 2 cm^^). 

The nature of the Ce2 bonding and the role of the 4/- 
states were also considered in CASSCF studies by Roos 
et al. in 0. As in it was found that the triple bond 
molecular configuration ((Tg)^(7ru)* couples in the 
molecular term. Two remaining nonbonding and local- 
ized 4/ electrons are in (/)-state (m = ±3). These two 
4/ electrons then can form the configurations (Afipg)'^, 
(4/(/)g)(4/0„), and (4/(/)„)^, which result in molecular 
terms with angular momenta = and 6. The six 
calculated low-lying levels again are almost degenerate 
having approximately the same spectroscopic parameters 
(without the spin-orbit interaction): = 2.63 — 2.66 
A, cva = 166 - 189 cm-i, = 2.61 - 2.68 eV (The 
CASPT2 level of theory and a basis set of VQZP qual- 
ity) . The spin-orbit coupling has only small effect on the 
computed properties, in particular it reduces to 2.51 
eV in good agreement with the experimental estimate 
Do = 2.57 eV [from the third law]. Calculations of Roos 
et al. give a longer Rg and a softer Wc in comparison 
with the values of Cao and Dolg [H. Both theoretical 
estimates of Wc 0j Q are smaller than the experimen- 
tal result [^: oje = 245.4 ± 4.2 cm^^ [the force constant 
ke = 2.28(1) mdyne/A]. 

Thus, the standard well-known many-electron proce- 
dures iSj] correctly describe the Ce2 chemical bonding 
but find six nearly degenerate low-lying molecular terms 
(iE+, lE", 3S-, 3E+, ^6g and ^6^) instead of a well de- 
fined ground state. The difficulty lies in the smallness of 
the splitting energies of these excitations requiring a spe- 
cial treatment. In the following we formulate a simplified 
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valence bond (VB) model accounting for the 4/ electrons 
in the cerium dimer, find the important interactions and 
obtain the "fine" energy spectrum of the electron exci- 
tations. Our approach is based on the picture of 4/-4/ 
interactions as rightly found in [4, 5j , but requires a more 
deep insight in the couplings which experience the 4/ 
electrons in the molecule. 



II. MODEL 

We start with the description of the triple bond of Ce2. 
The ground state electronic configuration of the cerium 
atom is (4/)i(5d)^(6s)^ with the ground term. This 
implies that the Ce atom is a rare case of the first Hund 
rule violation Q. The nature of such behavior lies in 
"quartic Fermi hole" which effectively separates the va- 
lence electrons and decreases their Coulomb repulsion ^ . 
Already at this stage of the investigation it is clear that 
the 4/ electron is not independent, through the Coulomb 
repulsion it is strongly coupled with the valence 5c? and 
6s electrons on the same cerium site. 

However, the (4/)^(5d)^(6s)^ configuration with the 
closed 6s— shell is not favorable for chemical bonding. 
Therefore the electronic atomic state of Ce is first trans- 
formed to the (4/)^(5(i)^(6s)^ configuration with the ^H^ 
term situated at 0.299 eV above the ^6*4 ground term 
0. In the (4/)^(5d)^(6s)^ configuration both 5d and 6s 
shells become open and form the triple bond with six va- 
lence electrons: (6scrg)^(5(i7r„)'*. The bonding state 
of 5d and 6s electrons (without the 4/ states) is well 
understood 0j in particular it occurs in La.2 where 
the molecular characteristics {Dq = 2.5 eV, fee = 2.28(1) 
mdyne/A [|, the (5d)i(6s)2 initial and (5ci)2(6s)i final 
[at 0.33 eV] configuration) are very close to those in Ce2. 
Although the triple bond is strong, the total energy gain 
is weakened by 0.6 eV because of the necessity to promote 
one 6s-electron to the 5(i-shell for both cerium atoms. 

We consider the first cerium atom at (0,0,0) and the 
second at (0,0, i?e) (so that the z-axis is the molecular 
axis) and find the valence density distribution of the four 
5d and two 6s valence electrons. The 5d—6s ground 
state implies an axially symmetric electron density. The 
conclusion can also be proved by considering the den- 
sity distributions of two s and four d-electrons in the 
(6s(Tg)^(5d7r„)'* configuration. 

Now to the 5d and 6s valence states we add two 4 
electrons. The main assumption supported by 0: 
is that to a first approximation the 4/ electrons do 
not contribute to the chemical bonding and remain lo- 
calized. Then it follows that each 4/ electron should 



be either in the Yi^'^{Q,ip) — C sin'^ 9 cos 3(/? or in the 



Yg '''(9,1^9) = Csin 9sin3iy9 angular orbital state where 
C — —y^70/n/8. (Here we deal with the real spherical 
harmonics y^™'* (t = c for cosine- and t = s for sine-like 
azimuthal (ip-) dependence) defined in Rcf. [l3|-) That 
is because in both these states the Coulomb intrasite re- 
pulsion between the 4/ electron and the bonding 5d and 




FIG. 1: The azimuthal dependence of two 4/ orbital func- 
tions: (a) cos-like Vj^g, (b) sin-like Vj^g. The z-axis which 
is perpendicular to the plane of the picture (not shown) is 
parallel to the Ce2 molecular axis. 



6s states at the cerium site is minimal (see Appendix A) . 
Indeed, in these m = ±3 /—states the angular density de- 
pendence is sin^ 9 with the maximum at 9 = 7r/2, while 
in the binding 5d states maximum lies at 9 = 7r/4 and 
in 6s states at 9 = (considering their overlap). Thus, 
the occupation of the m = ±3 /—states ensures a cer- 
tain space separation of the populated 4/ states from the 
binding 5d and 6s electrons, which effectively minimizes 
the intrasite 4/ — 6s and 4/ — 6s Coulomb repulsions. 
The situation is discussed in detail in Appendix A, no- 
tice that the conclusion about the occupancy of the 4/0 
states is also explicitly supported by the calculations of 
Roos et al. Q and implicitly by Cao and Dolg y. 

Thus, at each cerium site the 4/ electron space is lim- 
ited to two TO = ±3 orbitals, see Fig. [TJ 



■0/,c(9,cp,r) 



(la) 
(lb) 



Here Rf{r) is the radial part of the 4/ states. 

Since there is very little overlap between 4/ functions 
belonging to diff'erent sites, the valence bond (VB or 
Heitler-London) [ll|, [l^] approach should give the best 
description to the 4/ electronic problem. Its basis set 
consists of the two electron 4/ states: 

*o„,Oi,(^l,Sl2; ^2, S2z) = V'/.Oa (^l)wis '0/,Oi,(^2)u2s, (2) 

where the indices Oa and Ob refer to two orbital sites at 
sites a and 6, i.e. Oa,oi, = c or s (orbital functions defined 
in Eq. ([Talb) and visualized in Fig. [T|), s^ is the spin 
projection number (±1/2), Ug is the corresponding spin 
function (a or /3). 

The two electron wave function 
*o„,Oi,(^i,si2; ^2,S22)=*o(^i, ^2) *s(siz; S2z) IS a 



product of the orbital part (* 
be labeled as 



1 — 4), which can 



|c,s), |s,c), |c,c), |s, s), (3a) 

and the spin function j- (j = 1 — 4), 

(a, a), (a,/3), (/3,a), (/?,/?). (3b) 

(Here we omit cerium site indices a and b, that is |c, s) = 
\ca, Sb) etc.) Therefore, in total there are 16 two electron 
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basis states t {t =1-16), which are then antisymmetrized 
in respect to the electron permutation. 

In atomic units the corresponding valence bond Hamil- 
tonian [111, [13 given by 



2 2 ri2 



(4a) 



where is the spin-orbit coupling considered in detail 
below (Sec HrC]) . while 



C/(rO = C/„(r) + C/b(f). 



(4b) 



Here Ua{h){^) < is an effective 4/ attraction to the 
nucleus a (or b) screened by all core electrons. 

It is well known that the evaluation of the matrix ele- 
ments of H^^ in the space of the functions ^ leads to 
the following equation: 



(5) 



Here Eq — £^[(6scrg)^(5d7r„)''] is the background molecu- 
lar energy, 6 is the Kronecker delta, while J and K are 
the direct and exchange matrix elements which we con- 
sider in Sec. Ill Al and Sec. Ill Bl respectively. 

Taking into account the orbital degeneracy of 4/ states 
the matrix elements of _ff (16x16) where H = J, K ot 
H"" axe conveniently written in the spin block form. 



H = 



/Hii Hi2 Hi3 Hii" 

H21 H22 H23 i?24 

H31 i?32 -^33 H34 

\H41 H42 H43 H44J 



(6) 



where the indices si, §2 refer to the spin part, Eq. (j3bp . 
As usual, H2\ = etc. Each of the blocks in Eq. © 
can be further represented as a 4x4 orbital matrix of the 
following structure: 



Ml ci di ei^ 

ci ai d2 62 

d\ d2 a2 C2 

Ve^ 63 C2 a2/ 



(7) 



Here the matrix elements are evaluated between the two 
/—orbital sets [11,12, Eq. (I3al) ] with the fixed spin in- 
dices [si,S2, Eq. pb[) ]. In our basis with real spherical 
harmonics, Oj, hi, Ci are real, while di and are imag- 
inary values describing the interaction with an external 
magnetic field or the spin-orbit coupling. The specific 
form of the matrix elements is considered below. 



A. 4/ 4/ direct Coulomb repulsion 

Here we consider the matrix of the direct Coulomb 
repulsion: 

Jt^M = y V'oi,a(fl)V'o3,a(?l)^(?l,n2) 

ipo2.b{'^'2)'ipoi.b{'^'2)dfidr2, (8) 



where V{'ri,r2) = l/|n. — ^2! and the orbital indices 
(01,02) belong to the composite index ti, while (03,04) 
to t2, Eq. @. 

From symmetry arguments it can be shown (Appendix 
B) that there are three types of the matrix elements, 

A = (c, S\V\C, S) = (Ca, Ca\V\sb, Sb) = (Sa, Sa|l^|cb, Cfe), 
B = {c,c\V\c,c) = {Ca,Ca\V\cb,Cb) = {Sa, Sa\V\sb, Sb), 
C = (C, S\V\S, C) = (Ca, Sall^lCb, Sb) = (Ca, Sa\V\sb, Cb), 

(9) 

while 

(C, C\V\C, S) = {Ca,Ca\V\cb, Sb) = (Sa, Sa| V|c6, Sb) = 0. 

Here we use both standard and so called chemist's no- 
tations [l6j . However, it turns out that not all param- 
eters in ^ are independent. Analyzing the azimuthal 
dependence of the integrands and "the azimuthal selec- 
tion rule" (Appendix B), one finds 



B~-A^2C. 



(10) 



In the H = J case in Eq. ([6|) only the diagonal blocks 
Jss are relevant (Jsi,s2 7^ — ^2) because the 

matrix elements ([S]) are independent of the spin indices. 
Furthermore, the Jss blocks are identical, i.e. J°^^ = Jss, 
Eq. The J°^^ matrix is given by Eq. ([7]) where ai = 
A, a2 ~ B, ci — C2 — C and di = d2 = ei = €2 = 0. 



B. 'if — 'if exchange interaction 

For the matrix element of the 4/ — 4/ exchange inter- 
action we obtain 

ry- ( QO.,b rro.,h . rra,h na^h 

^t-i, t2 V 01, 04 ^02, 03 01, 04 ^02, 03 

'^01,04:02,03 j "si^, S4^ "S2i, S3^, V 7 

where the composite indices ti and t2 refer to 
(oi,a,siz; 02,6, S2^) and (03^^,53^; 04,^,342) as described 
earlier, Eq. ([2]), and 



dri>oiA^'^02,b{.^^ (12a) 



C^o;'o2 = j rfr^oi.a(r)C/(f)Vo2,fc(r). (12b) 

Finally, the Coulomb exchange integral in (jlip is given 
by [compare with Jti,t2 in Eq. ([5])] 



^ Oi,04;02,03 



'0oi,a(fl)'0O4,b(^l) V{ri,r2) 

i^02,bi^2)i>03.ai'^2)dridr2- (12c) 



If two first terms on the right hand side of pTjl differ from 
zero then they are leading, while the third (Coulomb) 
term is just a minor correction. 
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In our case there are two 4/ functions at each site, 
Eq. P^ - (|lb[) . but from symmetry it follows that in all 
overlap integrals S^i^^^^ there is only one nonzero value, 

S = {Ca\cb} = {Sa\sb) > 0, (13) 

while {ca\sb) = 0. Analogously, Uexc = {ca\Ua\cb) = 
{sa\Ua\sb) and {ca\Ua\sb) — 0. It is convenient to intro- 
duce the notation 

SE^-2SUexc>0 (14) 

for the leading term. 

For the exchange Coulomb interaction V'^^'^, Eq. (|12cp . 
one finds three types of integrals with nonzero values 
(compare with Eqs. 

A' = (c, S\V\C, S)exc = {Ca,Cb\V\Sa, Sfc), (15a) 
B' = (c, c\V\c, C)exc = {Ca,Cb\V\Ca, Cfc), (15b) 
C = (c, s\V\s, C)exc = {Ca,Sb\V\Ca, Sfc). (15c) 



TABLE I: Calculated 4/ direct {A,B,C) and exchange 
(A', B' , C') Coulomb integrals, the overlap integral S, the VB 
attraction Uexc and the 4/ exchange parameter SE for vari- 
ous Ce-Ce distances R (in A). All energy parameters are in 
Kelvin (K). 

" Ce2 equilibrium distance for ^Eg , in Ref. [j, 
Ce2 equilibrium distance for ^E+, ^6g, '^E^, ^6ti in Ref. [^, 
Ce2 equilibrium distance for ^E~, '^E~ in Ref. [3]. 

R =3.43 A and 3.65 A correspond to the nearest neighbor 

distance in a-Ce and 7-Ce [l^, respectively. 





R =2.62" 


2.63'' 


2.66'= 


2.70 


3.43 (a) 


3.65 (7) 


B 


2324.45 


2286.28 


2174.68 


2032.49 


643.73 


471.58 


A 


2323.98 


2285.82 


2174.26 


2032.12 


643.69 


471.56 


C 


0.24 


0.23 


0.21 


0.19 


0.02 


0.01 


B' 


1.82 


1.75 


1.56 


1.34 


0.09 


0.04 


A' 


1.67 


1.61 


1.43 


1.23 


0.09 


0.04 


C 


0.07 


0.07 


0.06 


0.05 








S X 10^ 


4.224 


4.148 


3.928 


3.654 


1.029 


0.714 


U exc 


-1727.81 


-1688.96 


-1578.02 


-1442.24 


-309.24 


-200.34 


5E 


-14.6 


-14.01 


-12.4 


-10.54 


-0.64 


-0.28 



Furthermore, there is the same relation as for the direct 
matrix elements A, B and C, Eq. ([TU|. see Appendix B: 



B' -A' = 2C", 



(16) 



Then, in the full 16x16 matrix H = K, Eq. ©, the 
exchange coupling is effective in the blocks Kl'^'^ (between 
states a, a), i^TIf^ and K^2'^ (between the states with a, 
[} and /?, a) and i^lf^ (between states /?, /?). Each of the 
blocks has the following (orbital) form: 



T/'exc 



/ -C" 5E-A' 

5E-A' -C 

SE-B' -C 

V -C SE^B'' 



.(17) 



(Here as in Ref. [ll| we have taken into account the effect 



of the half canceling of the 6E-terms in the final expres- 
sion due to the nonorthogonality of the /-functions.) 



C. Spin-orbit coupling 

Now in the framework of our model we consider the 4/ 
spin-orbit interaction: 



(Lisi +CL2S2, 



(18) 



where C = C/ ~ 862 K. This value corresponds to the 
Aso = (7/2)C =0.26 eV splitting between j = 5/2 and 
7/2 one electron atomic 4/ states. (The energies of the 
states were obtained by solving the spherically symmetric 
Dirac equation for atomic cerium states in the framework 
of the self-consistent field DFT-LDA atomic model.) In 
our model the 4/ quantum space is reduced to two orbital 
states with niz ~ ±3 at each site, and H''° transforms to 



(19) 



Therefore, the nonzero spin blocks of the full matrix are 
^s°8 — 1 ~ 4), Eq. (|n]), with the following general 
structure of each block (see Appendix C): 



nii 7i2i' 

mil 

—nii —mil 

-71,2* —m2i 



(20) 



-1, in the 



in the 

= — 1 and in block 



In the Hl° block ni — mi — 1, 712 ^ m2 - 
block ni = n2 ~ —1, mi = 1112 = 1 

block Til =712 = 1, "nil = 7712 
77,1 = 7771 = —1) n2 = 7772 = 1- 

The spin-orbit coupling, Eq. splits the four single 
electron states at each cerium site in two doublets: the 
lowest level with = ±5/2 and the energy i?(±5/2) = 
— 3C/2 has the eigenfunctions 7/'5/2,+ ~ exp(377iz(p)/3, 
"05/2,- ~ exp(— 377iz(/?)a. The other level with jz = ±7/2 
and the energy i?(±7/2) = ±3C/2 has the functions 
"07/2,+ ^ cxp(3777,7(p)Q!, 'ip7/2,~- ^ exp(— 377i7</?)/3. Cor- 
respondingly, two 4/ electrons have the following energy 
spectrum: Ego.i = — 3C (4 degenerate states), Eso,2 = 
(8 states), Ego.s = 3C (4 states). The inclusion of other 
interactions modifies the picture. 



III. RESULTS AND DISCUSSION 

All interaction parameters of the model introduced in 
Sec. |TT]were evaluated numerically with our valence bond 
program p^ . Table HI There, we have used the overlap- 
ping cerium relativistic atomic 4/ functions (the large 
j = 5/2 component for Rf{r)) calculated within the den- 
sity functional theory (DFT) in the local density approx- 
imation (LDA) 

Below we discuss the results of the valence bond model 
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with (i?™ ^ 0) and without {H^° = 0) the spin-orbit 
couphng. 

A. VB without the spin-orbit coupling 

From Table |T] we see that the largest parameters are 
A and B which correspond to the direct Coulomb repul- 
sion, Sec. Ill Al Therefore, as a first approximation we 
consider only the matrix of the direct Coulomb repulsion 
J. In that case from the structure of the spin matrix it 
follows that the triplet and singlet states coincide. Hav- 
ing diagonalized the corresponding 4x4 orbital matrix 
analytically, we find its eigenvalues, 

Ex- A = -C, (21a) 
Ei,z -A = C (21b) 
Ei- A = 3C. (21c) 

(The eigenvectors are quoted in Eqs. ()C7a|) - (|C7cp . Ap- 
pendix C.) We conclude that the energy splittings are 
determined by the parameter C, Eq. ([U]), which is very 
small, while A should be incorporated into the back- 
ground molecular energy Eq = E[{Qs(jgY{bdTTy)'^]. 

Notice that the eigenvalues E2 and E^ are degenerate 
which implies a nonzero orbital momentum. In Appendix 
C we show that indeed the two components are attributed 
to Mz = ±6. Taking into account other properties of the 
eigenvectors one can show that the molecular terms are 
identified as S^, 6g and S^, respectively. However, the 
presented picture is not full because the integral C is very 
small. Table H] If C were of the order of A, the direct 
Coulomb splitting would have been dominant, but since 
C ^ 0.2 K we should consider other relevant interactions. 

Now in addition to the direct Colomb repulsion J we 
consider the 4/ — 4/ exchange interaction K, Sec. Ill Bl A 
small overlap S between two relevant 4/ functions splits 
the singlet and triplet states. If there were only one or- 
bital for each site, the spin singlet state would have been 
lower in energy than the spin triplet state. Indeed, in the 
single orbital valence bond problem [ll|, [l^] we have 

Es=Eo- 6E, ET = Ea+ SE, (22) 

where 6E is given by Eq. ((Ti| . In our case of two 4/ or- 
bitals at each site the reverse is true. This can be clearly 
seen if for the time being we ignore the exchange repul- 
sion {A' = B' = C — 0), and treat the electron site 
exchange as a perturbation. As an example we consider 
the state (j21al) with the eigenstate (|C7ap . Its orbital part 
can be written as [~Ca{ria)sbir2b) + Sa{ria)cb{r2b)]/V2, 
where fia = \ri - Ra\, rib = \r2 - Rb\- In the spin 
singlet state we then symmetrize the orbital part, ob- 
tain [-Ca{ria)Sb{f2b) - Ca{r2a)sb{rib) + Sa{ria)cb{r2b) + 

Sa{r2a)cb{fib)]/2 and find E{^J:-) ^ Ei + SE. In the 
triplet state we antisymmetrize the orbital part, ob- 
tain [-Ca{ria)Sb{r2b) + Ca{f2a)Sb{rib) + Sa(ria)Cb{r2b) " 

Sa{r2a)cb{rib)]/'2 and find ^(^S;) = E^ - SE. Thus, 



TABLE II: The VB energy spectrum of Ce2 without the 
4/ spin-orbit coupling ((" = 0) and the exchange Coulomb 
repulsion (A' = B' = C" = 0). 



singlets 


E-Eo 


triplets 


E-Eo 




-C + SE 




-C -5E 




C ~5E 




C + SE 




3C-SE 




3C + 5E 



TABLE III: The analytical expressions for the energy levels 
of Ce2 in the VB model with the spin-orbit coupling {A' — 0, 
B' = 0, C" = 0), P is the state parity. 





E 


deg. 


n(J.) 


A(L.) 


E(5'.) 


P 


1 


~^y9C +4C'^ + C -3E 


1 











9 


2 


-^/QC+AC^ + C + SE 


1 











u 


3 


-3C + C + 5E 


2 


T5 


T6 


±1 


u 


4 


-C-5E 


2 


Tl 





Tl 


9 


5 


C-5E 


2 


T6 


=F6 





9 


6 


C + SE 


2 


T6 


=F6 





u 


7 


iC + SE 


2 


Tl 





Tl 


u 


8 


yjQC^ + AC^ + C ~ SE 


1 











9 


9 


3C + C + SE 


2 


+7 


T6 


Tl 


u 


10 


+ ^C"^ + C + SE 


1 











u 



the Ce2 ground state is the triplet state in contrast 
to the standard (diatomic molecule) case. The origin of 
this anomalous behavior is most likely connected with the 
fact that the orbital part of the triplet state is even. The 
results for all states are summarized in Table HIl which 
holds even if SE is not a perturbation. The energy levels 
are also schematically shown in Fig. [51 

Thus, the situation in Ce2 could be an interesting ex- 
ample of the triplet ground state. In practice however 
the relatively strong (in comparison with SE) spin-orbit 
coupling changes the electron spectrum. 

B. VB with the spin-orbit coupling 

In this section we discuss the most general case, when 
the full valence bond Hamiltonian is given by Eq. 

(Sil). 

Surprisingly, an analytical solution can be figured out 
even for this case except for the small Coulomb repulsion 
which accompanies the 4/ exchange transitions [i.e. A' — 
0, B' = 0, C" = 0, Eq. (fT5a|) - ((T5c)) ]. The analytical 
expressions for the energy levels are quoted in Table IIIII 
Notice that the three lowest states are those found in Ref. 
4] [{VI = 0)g, 0„ and 5„] and Ref. Q (Table 6). 

For completeness we can include the 4/ electron ex- 
change repulsion {A' , B' , C) and compute the energy 
levels numerically. The calculated energy values for char- 
acteristic internuclear distances R and the parameters of 
Tabled] are quoted in Table |lVl Since A', B' and C" are 
very small, they introduce only a minor correction to the 
analytical expressions. 
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TABLE IV: Calculated energy levels (in K) of Ce2 at various 
characteristic internuclear distances R (in A) with the spin- 
orbit coupling, see also notes for Table H] 





R =2.62 


2.63 


2.66 


2.70 


3.43 (a) 


3 65 ('t'I 


X 


-9602 1 1 


-2601 46 


-2599 69 


-2597 64 


-2586 71 


-2586 31 


2 


-2'i6Q 42 


-2570 ns 

ZjO I U . yJO 


-2571 90 


-2573 99 


-2585 25 


-2585 67 


3 


-2569.42 


-2570.08 


-2571.90 


-2573.99 


-2585.25 


-2585.67 


4 


-16.59 


-15.92 


-14.11 


-12.02 


-0.75 


-0.33 


5a 


-16.26 


-15.60 


-13.82 


-11.75 


-0.71 


-0.31 


5b 


-15.96 


-15.32 


-13.56 


-11.53 


-0.71 


-0.31 


6a 


16.44 


15.78 


13.98 


11.91 


0.75 


0.33 


66 


16.74 


16.06 


14.24 


12.13 


0.75 


0.33 


7 


17.07 


16.38 


14.53 


12.40 


0.79 


0.35 


8 


2569.90 


2570.54 


2572.32 


2574.37 


2585.29 


2585.69 


9 


2602.59 


2601.92 


2600.11 


2598.02 


2586.75 


2586.33 


10 


2602.59 


2601.92 


2600.11 


2598.02 


2586.75 


2586.33 



Usually, for labeling of molecular states one uses the 
standard many-electron spectroscopic notation ^'^"'"-'^Ao 
[l7| . However, in the — Sz = states the inclusion 
of the spin-orbit coupling mixes up singlet (S* = 0) with 
triplet {S = 1) states and the partial weights of singlet 
and triplet contributions in the final Sz = states are 
nearly equal. (The effect was found in 0, In our 
model for the lowest three states we obtain: 

^-1 = a*(3l];; Sz=0)- 6*(iS+; Sz - 0), (23a) 

*2 = -a*('S-; Sz=0)+ S+; Sz = 0), (23b) 

^-3,1 = -a*('s;; Sz = l) + Sj; Sz - 1), (23c) 

^-3,2 = b^i^i:-; Sz = -1) + Sj; Sz = -1). (23d) 

Here a = 0.7072, b = 0.7070, ^(^Sg ; Sz = 0) is the 
Sz — component of the "^S^ triplet etc. Notice that the 
and states in vj/g form the level (f2 = 5), 
see FigH] below. In the following we will use the short 
notation (1,3) for the spin multiplicity (25* + 1) of the 
strongly mixed states. For example, the ground state is 
labeled as ^'^So,g (or ^'^So,g) etc. 

The exchange Coulomb interaction {A',B',C' ^ 0) 
splits the doublets 5 and 6 (fi = A = 6, S = 0, g,u), 
Table imi 

^■'66(,,„)^i'3l]+^^„)+i'3S-(^^„). (24) 

In the splitted lines (5a, 5b, 6a, 6b in Table IIV[) = 0, 
A = and E = 0. Although the ordering of three lowest 
levels is the same as in Ref. Q , the energy splittings are 
different. At R = 2.62 A we have E2 ~ Ei = E3 - Ei = 
32.7 K, while in the spin-orbit calculation of Cao and 
Dolg |4j E2- Ei<l cm-i and E3 ~ Ei < 2 cm^^. 

It is worth noting that the levels 2, 3 and 9, 10 are 
virtually degenerate. According to Table IIIII the energy 
difference between them is 2C^/3C < (C, SE). 

If the spin-orbit coupling is dominant (which is the 
case for Ce2), the twelve energy levels are clearly divided 
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FIG. 2: The schematic energy splittings in the VB model in 
the absence of the spin-orbit coupling: (A) only SE 7^ 0; (B) 
SE, C / 0, A',B',C' = ; (C) SE, A',B',C' / 0, C = 0; 
(B-f C) SE, A, B, C, A', B', C' / 0. In the case (C) the levels 
^E~, ^E^ and '^Ej, ^E~ are "accidentally" degenerate. 

in three groups: 

(1) ''^Eo(g), i'3So(„), '65(„). (^ E,o,i) (25a) 

(2) ^Ei(g^_), ^■3l]o(g,+), ^'^Eo((;,-)i 

(3) l'35^o(g), '67(„), I'^Sp^^j. (^ ^.0,3) (25c) 

Here the levels are given in the increasing energy order as 
in Tables Hill and HVl with some additional characteristics 
(the parity and the mirror reflection, when possible) in 
parentheses. If the other interactions are omitted (i.e. 
SE = 0, C = 0, and A' = B' = C = 0) then the energy 
spectrum is reduced to only three levels {Ego.i = — 3C, 
Eso,2 = 0, Eso,3 = +3C) discussed in Sec. Ill CI 

The relation between three important characteristic 
parameters of Ce2 is the following ( ^ SE ^ C. How- 
ever, the analytical and numerical solution can be used 
to study the model for any parameter set. It is clear that 
such a scenario is not relevant for Ce2, but it can pos- 
sibly occur in other electronic systems with localized / 
electrons. In particular, it is instructive to consider the 
case of a small spin-orbit coupling: (C, SE) ^ For 
that relation (provided that A' = B' = C = 0) one can 
still use the analytical solution. Table iHll but of course 
the order and level grouping will be different. The split- 
tings and relations between energy levels in the absence of 
the spin-orbit coupling (C = 0) are schematically shown 
in Fig. H 

In Fig.[2]the most realistic subcase (all relevant interac- 
tions included) is (B-l-C). If on top of it a small spin-orbit 
coupling is introduced, then the spin triplet states there 
are split according to the following scheme: 

^^Eo (<,,-)+ (^^_), (26a) 

^E- -^3Eo(„^_) -1-365 („), (26b) 

^''j:oiu,+)+%iu), (26c) 

2'j:+ ^^^oiu,+)+'^iiu,+)- (26d) 
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IV. CONCLUSIONS 



For the cerium dimer we proposed a model of two inter- 
acting 4/ electrons which can explain and refine the low- 
lying energy spectrum of Ce2 observed in many electron 
computational studies [3, § ■ The model is a modified va- 
lence bond approach where two relevant 4/ orbital states 
at each cerium site are explicitly taken into account, Eqs. 
([Ta)) . (|Tb)) . and Fig.[T] The other valence 6s and 5d states 
are grouped in a conventional triple chemical bond with 
the {6sag)'^{5dTTu)'^ closed molecular shell. 

In the framework of the 4/ model we can clearly trace 
the origin of the energy levels and establish intrinsic rela- 
tions between them. For the most important interactions 
the problem of finding the 4/ energy spectrum is solved 
analytically. Apart from the spin-orbit coupling =862 
K) the largest interaction is the direct Coulomb repulsion 
of the two localized 4/ electrons (described by parameters 
A, B) but the effective Coulomb parameter responsible 
for the energy splitting is very small (C ^0.2 K), Table 
m The next important interaction is the valence bond 
exchange of the two localized 4/ electrons {5E ~14 K) 
which is also accompanied by small exchange Coulomb 
repulsion {A',B' -1.5 K, C -0.07 K). 

As a result, the even ground state level with = 
A = E = (composed of ^^o(g) ^nd ^^Q(g] states, 
Eq. (I23ap ') is separated from first excited levels with 
O = A = E = (composed of '^So(„) and ^So(«) states, 
Eq. (|23bl) ) and ^Q^^u) by an energy gap of ~30 K. The 
ordering of these levels is in agreement with calculations 
of Cao and Dolg although the energy splittings are 
different. The Jz — ^ — ±5 excited states support mag- 
netic moments, giving rise to the magnetic susceptibility 
dependence x shown in Fig. [3] The magnetic suscep- 
tibility clearly displays two different regimes. At large 
temperatures (T > 30 K) it follows the Curie law indica- 
tive of "free" magnetic moments. At small temperatures 
(T < 20 K) X deviates from the Curie law and rapidly 
decreases as T — 0. Such low T behavior and disap- 
pearance of the "free" 4/ magnetic moments of cerium 
is caused by the population of the nonmagnetic ground 
state (n = A = E = 0) at the expense of the excited 
magnetic ones. The conclusion can be checked by exper- 
imental measurements. 

We have also studied the energy spectrum in the ab- 
sence of the spin-orbit coupling. Fig. [21 Here, the unusual 
feature is the ^Tj~ triplet ground state. 

The peculiarities of the electronic structure of the sim- 
plest cerium bond possibly can clarify the picture of 
phase transitions in solid cerium [Tsj . 

I am grateful to K.H. Michel, B.N. Plakhutin, A.V. 
Tsvyashchenko, E.V. Tkalya, A.V. Avdeenkov and B. 
Verberck for valuable discussions. 




TEMPERATURE (K) 

FIG. 3: The calculated magnetic susceptibility of the cerium 
dimer. The magnetic field is parallel to the molecular axis. 
The Curie law tail at high temperatures can be fit with the 
effective Bohr magneton number 2.42 /xb (about 300 K) per 
cerium site. 

Appendix A: Appendix A 

According to the recent many electron calculations 
(CASSCF) of Ce2 [3] each cerium site has a single 4/ 
electron which is virtually nonbonding and preserves its 
localized nature. Furthermore, without the spin-orbit 
coupling it occupies one of the two orbitally degenerate 
4/ states with the highest azimuthal angular momen- 
tum (m = ±3). Although it is not mentioned explicitly, 
the same orbital degeneracy follows from the low-lying 
molecular term spectrum computed by Cao and Dolg Q . 
Below we show that the lowering of the m = ±3 states is 
caused by the Coulomb repulsion between the 4/ states 
and the bonding electron states. 

We consider a simplified model with two bd and one 
6s bonding electrons at each cerium site and use the 
method of multipole interaction [l^. According to 
[1, Q one 6s electron and two hd electrons from each 
Ce site contribute to the (6scrg)^(5d7r„)^ triple bond. 
That implies the occupation of the y^^2 ^2" '3''" 
bital d-states (m = ±1). The resultant electron density 
{{Y^'^iQ, ip)f + {Y^''{<d, ip)f) is axially symmetric. Ex- 
panding it in the multipolar series one finds that the first 
nonspherical component responsible for splittings is given 
by the Y^^2 angular dependence of the Sd-state density, 

Pd,i=2{r. e, ^) = 2C,V2 Rlir) Yl,{Q, ^), (Al) 

where i?d(r) is the 5d atomic function radial dependence, 
and C2 = 0.090112 is the angular density matrix ele- 
ment [l9| , 

Cli^2 = (2, ls|2, 0|2, Is) = (2, lc|2, 0|2, Ic) 

= / dnYi^2imi=2imi=2m- (A2) 

Here f2 stands for the polar angles {Q, (p). Now one can 
calculate the matrix of the Coulomb interaction of the 
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TABLE V: Coefficients C° ^^2(1") c) or (m, s)] and cal- 

culated splittings i54/(r) of the 4/ states (in eV) caused by the 
quadupolar Coulomb repulsion with two occupied 5d states 
{md = ±1), see text for details. 

r (m = 0) (l,c), (l,s) (2,c), (2,s) (3,c), (3,s) 

Cj_f=2(T) 0.16821 0.12616 -0.21026 

£;4/(r) 0.21540 0.16156 -0.26925 



5d density with the 4/ electron. (In principle it will have 
the £ = 2 and £ ^ 4 contributions but the ^ = 4 
term is small and in the following will be ignored.) The 
matrix is diagonal in the space of the 4/ orbitals with 
the diagonal elements given by 

Vf^r = VdfC%^2{r)- (A3) 

Here r = (m, c) or (to, s) (m = — 3) is the azimuthal 
dependence of the 4/ orbital functions, Cji^2{''') = 
(3,r|2,0|3,T) [compare with ^M^], while 

%=2C,V2«2,<i/, (A4a) 
^^Mf^^J Jdrdr'Rlirf^R}{r'). (A4b) 

r< (r>) is the largest (the smallest) value of r and r' . 
V2,df describes the quadrupolar {£ = 2) component of 
the Sd — 4/ intrasite Coulomb repulsion [l^l- With 
the 5d and 4/ atomic radial dependencies (DFT-LDA 
atomic calculations 14]) we obtain V2,df — 7.1054 eV 
and Vdf = 1.2806 eV. The C^^^2('^) coefficients and 
the corresponding 4/ energy splittings are quoted in 
Table |Vl (Notice that due to the z— axial symmetry 
C/,f=2('™'C) = Cl^^2{m,s).) 

Thus, we conclude that the hd — Af Coulomb repulsion 
favors the lowering of the to — ±3 /—states. However, 
there are two other interactions which act in the opposite 
direction: the screened attraction to the other nucleus 
and the bd — 4/ exchange interaction. 

The other nucleus attraction leads to an axially sym- 
metric Coulomb potential which has the form of (|A3p . 
but with an interaction constant of the opposite sign. 
The bd — 4/ exchange interaction is more complex. In 
terms of the intrasite Coulomb multipole repulsion it cor- 
responds to the odd (€ = 1,3,5) transition density con- 
tributions, and the dipolar {t — 1) contribution is rather 
strong [19]. However, this exchange interaction leads to 
the high spin {S — 2) ground state (the ^H^ term of the 
Afbd^ds configuration) of the cerium atom, which in the 
real molecule is effectively suppressed by the triple bond 
formation with a zero spin state of the binding 5d and 6s 
states. Therefore, the bond formation mechanism dimin- 
ishes the 4/ — 5d exchange interaction, while the 4/ — 5d 
Coulomb repulsion (IA3|) remains and lowers the to = ±3 
/—states in respect to the others [1,[H]- 



Appendix B: Appendix B 

The relevant Coulomb integrals, Eq. ([9]), can 
be simplified by considering their azimuthal parts, 
i.e. the azimuthal dependence of the corresponding 
ij^oi.airi) ipo3,a{ri) functions standing on the right and 
their counterparts V'o2,&(^2) 4'o4,,b{r2) on the left from the 
Coulomb interaction V{'ri,'r2) = 1/'^12, Eq. In the 
following we define = ri — Ra and ri, — r2 — Rb- The 
important observation is that the Coulomb integral dif- 
fers from zero only if the azimuthal density dependencies 
on the left and on the right side have the same functional 
form. This conclusion can be derived by various ways, in 
particular it follows from the multipole expansion of the 
Coulomb interaction 

V{Ta.n) = Y.Vl{Ta.n)Yafa)Yah). (Bl) 

Here the angular functions Yj^{fa) and Y^{fi,) have the 
same azimuthal dependence (r) which is a trivial conclu- 
sion if they refer to the same site: Ra ~ Rb- In case of 
two different sites (aligned along the z— axis) one can ex- 
pand a spherical harmonic centered at one site in terms 
of the spherical harmonics centered at the other site fl^ . 
It turns out that the procedure conserves the azimuthal 
dependence of the spherical harmonics translated along 
the z— axis ^13*], which proves our initial statement. 

In the following in deriving relations between the 
Coulomb integrals, Eq. we use chemist's notations 
for the matrix elements 16]. Further, we explicitly sin- 
gle out the azimuthal dependence , 03 (</') out of the 
density functions il}oua{ra)'4>o:i.a{'^a) with oi, 03 = c, s, 
Eq. ([Talb). We then arrive at three different dependen- 
cies: 

PcAv) = {cosiff = ^ + ^ cos6(p, (B2a) 
Ps,3((p) = (sin3.^)^ = i - icos6<^, (B2b) 

Pc,s{f) ^ cos3(y9 sin 393 = — sin6iy9. (B2c) 

From the azimuthal selection rule it follows that 
{ca,Ca\V\cb,Sb) = £> = 0, Eq. (H]) and (HU]), and 

{Ca,Ca\V\cb, Cb) = (s^, Sa|y|s6, Sfc) = B. 

In order to prove the relation (jlOp we consider the 
structure of the general Coulomb integral, 

{la,3a\V\2bAb) f dnf{e)lA{Ma\Mb), (B3) 

where Q stands for {Qa,Qb,ra,rb), while /^(MolMb) is 
the azimuthal integral with two identical (/^-dependencies 
Ma and Mb. For example, lA{Qa\Qb) with Ma = Mb = 
M — 6 stand for the azimuthal integral between cos6(/5 
and cosGip or between sin6(/3 and sin 695: 

lA,c{Ma\Mb) = / dLpa I difib 

Jo Jo 

cos{Mipa)V{ra,rb)cos{Mipb). (B4) 
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Each term on the right hand side of (|B3I) has the following 
form: 

f dnf{e)lA{Ma\Mb)= [sinQadea I rldra 
Jn J J 

j sinSbdeb J rldn C [P|(cos ea)P|(cos eb)]^ 

R){ra)R]{rb)lA{Ma\Mb). (B5) 

Here PJ^^(x) is the associated Legendre polynomial 
[To| , C = — -^/TO/tt/S, Rj{r) is the radial dependence 
of the 4/ functions. Taking into account three relevant 
azimuthal dependencies (|B2a|) - (IB2c[) we get 

{Ca,Ca\VWsb) = \ f dnf{n)ilAiOa,Ob) ~ lAi6a\6b)), 
{Ca,Ca\V\cb,Cb) = ] / ^^(/(r!) (/^ (Oa , Ofc) + (6a 1 6b)) , 

4 Jn 

iCa,Sa\V\cb,Sb) ^ ] [ ^^(/(r!) /a (6a 1 6b) . (B6) 

4 Jn 

From Eqs. (jB6|) we immediately obtain 

(Ca,CQ|y|Cfc,C6) - (Ca,Ca|y|s6,Sfc) = 2 (Ca , Sa | ^| Cfe, Sfc) , 

(B7) 

or B- A==2C, Eq. 

For the exchange charge Pabir) = il^ouair) i^o^.bir) 
(oi, 02 — c, s) which corresponds to the 4/ electron tran- 
sition from the site a to the site b we again have the three 
azimuthal dependencies, Eqs. (jB2aP - (IB2cp . Therefore, 
proceeding as before albeit with a modified expressions 
for the integrals we obtain the analogous relation 

(Ca, Cb\V\Ca,Cb) - (Ca, C6|y|Sa, Sfc) = 2(Ca, Sh|y|Ca, Sfc), 

(B8) 

or B' -A' = 2C', Eq. (HH). As before, the result is based 
solely on the azimuthal selection rule. 

Appendix C: Appendix C 



where /is is the Bohr magneton. 

Further, the 16x16 spin matrix Sz^ Eq. has two 
nonzero diagonal blocks 5*2(1, 1) and 5^(4,4) of the form 

5,(s,s) = ±Msi, (C3) 

where the plus sign refers to the a, a block (s, s = 1,1) 
and the minus sign to the /3, /3 block (s,s = 4,4), while 
1 stands for the unit 4x4 matrix. 

The molecular interaction with an external magnetic 
field H is given by 

V^ag = -MH, (C4) 

with the molecular magnetic moment 

M^2S + L. (C5) 

For the magnetic field H applied along the z— axis we 
obtain 

y^as = -M25. + L.), (C6) 

where h = HfiB- The relations ([02]), jCSj and (IC6| are 
used to calculate the z-component of the orbital, spin 
and full momentum for each two-electron level. 

For example, by means of Eq. (jC2[) we can find the 
orbital quantum number required for the characterization 
of the degenerate states, which are considered in Sec. 
Ill Al There we have discussed the three energy levels 
given by Eq. (|21ap - (l21cp . The corresponding eigenvectors 
are 

(-l,l,0,0)/\/2, (C7a) 
(1,1,0,0)/V2, (0,0,-l,l)/V2, (C7b) 
(0,0,1,1)/n/2. (C7c) 

In the presence of the external magnetic field the 
Hamiltonian V{f — f) of the direct Coulomb repulsion 
(see Sec. Ill A|) is modified: 

H^Ebond + Vif-f) + HL, (C8) 



The action of the z— component of the orbital momen- 
tum on the real spherical harmonics 1^™''^, yj*"'" [lOj (the 
z-axis coincides with the molecular axis) is given by 



(Cl) 



The matrix elements of — Li,z + L2,z in the space 
of two electron 4/ states (Eqs. ([3a|) and ([3b|) can be 
obtained from these relations. Explicitly, the matrix of 
Lz has the nonzero blocks Lz(s, s) (s = 1 — 4) [see Eq. 
®] of the form 



Lzis,s) = fJ.B 



( 


V Zi 



3i -iis 

ii -ii 

-■ii 

3i 



(C2) 



where the last term is the orbital polarization given by 
i/i = —Lzh (the first term Ebond is the binding energy 
of Ce2 which is not important for our discussion) . From 
(jCSP we obtain four eigenvalues (i = 1 — 4) 

E, -Eo= -C, C - 6h, C + 6h, iC, (C9) 

and the corresponding eigenvectors: (— 1, 1, 0, 0)/\/2, 
(i,i,-l,l)/V^, (i,i,l,-l)/%/2, and (0, 0, 1, 1)/V^. 

Notice, that at zero magnetic field {h — 0) i?2 = E3 
as in (j21b|) . On the other hand, from (|C9|) it is clear that 
two components correspond to Mz = ±6 and therefore in 
this doublet Lz = 6. In other two non-degenerate states 
Lz — Mz = (i.e. they are E-states). 
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